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We study the propagation of neutrinos in gravitational fields using wave functions that are exact 
to first order in the metric deviation. For illustrative purposes, the geometrical background is 
represented by the Lense-Thirring metric. We derive explicit expressions for neutrino deflection, 
helicity transitions, flavor oscillations and oscillation Hamiltonian. 
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I. INTRODUCTION 

Because of their coexistence in a host of problems, the interaction of neutrinos with gravitational Helds is of 
interest in astrophysics and cosmology. At the same time, neutrino detectors of increasing sensitivity and scope 
require a better understanding of propagation and oscillation properties of neutrinos in the neighborhood of massive 
astrophysical objects. 

While in a few instances the treatment of gravity may necessitate the complete apparatus of general relativity, using 
a weak field approximation is sufficient in most problems that deal with neutrinos. Even in the case of black holes, a 
first order approximation in the metric deviation may still be adequate up to distances of a few Schwarzscliild radii 
from the gravitational source. In the majority of instances, gravitation is treated as a quasi-classical external field. 
Even so, the calculations tend to be complex. 

Useful, simplified treatments of the problem of neutrinos in gravitational fields can be found in the literature, 
notably in the heuristic approach of Cardall and Fuller [1] in which the authors construct an effective mechanical 
four-momentum operator that incorporates spin and matter effects. 

A more direct approach would require solving the Dirac equation completely which is in general very difficult. It was 
however shown in [2,3] that the covariant Dirac equation can be solved exactly to first order in the metric deviation. 
This procedure is manifestly covariant and gauge-invariant and only requires the evaluation of path integrals that are 
trivial for most physically relevant metrics. It can also be extended to include electromagnetic fields and the effect of 
media on the propagation of the particles. These are all issues of principle one meets in the interaction of quantum 
systems with semiclassical gravitational fields [4]. They have been raised time and again in the literature [5]. We think 
it advantageous to address them within the context of an approach that is consistent and reproduces all gravitational 
effects that are known, or are supposed to apply to fermions [6,3,7]. We propose to apply this method to the study 
of neutrino optics, helicity transitions and flavor oscillations. 

In Sec. II we briefly review the essential points of the solution found in [2], The neutrino geometrical optics is given 
in Sec. Ill and applied in particular to the propagation of neutrinos in a Lense-Thirring background for the distinct 
cases of propagation parallel and orthogonal to the angular velocity of the gravitational source. Helicity transitions 
and flavor oscillations in the absence of matter are dealt with in Sec. IV and Sec. V for the same two directions of 
propagation. Sec. VI contains the conclusions. 
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II. THE COVARIANT DIRAC EQUATION 


The behavior of spin-1/2 particles in the presence of a gravitational field g^ is determined by the covariant Dirac 
equation 

[i^ixyDn - m]^{x) = 0, (II. 1) 

where V M is the covariant derivative, T M (a:) the spin connection and the matrices 7 P ( x ) satisfy 

the relations {^(x), 7 ^(a:)} = 2 g^ v . Both T Al (a;) and 7 ^( 2 ;) can be obtained from the usual constant Dirac matrices 
by using the vierbein fields e/ and the relations 

7^0) = e &( x )l & , r M (a;) = , (II.2) 

where <j a0 = f [ 7 °/7^]- A semicolon and a comma are frequently used as alternative ways to indicate covariant and 
partial derivatives respectively. We use units H = c = 1 throughout the paper. 

Eq. (II.1) can be solved exactly to first order in the metric deviation j Ml/ (x) = g— rj^, where the Minkowski 
metric rj^ has signature -2. This is achieved by first transforming (II.1) into the equation 

[« 7 I/ (a;)V !/ — m]^{x) = 0 , (II.3) 

where 

$(*) = S’ -1 * (s), S(x) = , $,(*) = V J dz x T x {z), ^(x) = S~ 1 ^(x)S . (II.4) 

By multiplying (II.3) on the left by (— 47 " (a;) — m), we obtain the equation 

GT'V^+m 2 )'^ ®)=0, (II-5) 

whose solution 

#(*) = (II. 6 ) 

is exact to first order. The operator 4 >g( 2 :) is defined as 

®G = -^J dz x [ 7 a x,p(z) -'y f 3x, a ( z )]^ al3 ( z ) + \j dz x j a xk a , (II.7) 

[L a0 (z), T 0 (x)] = ((*“ - z a )kP - (x 0 - z 0 )k a ) Mx), [k a , tfo(z)] = id a *o , 

and To (a;) satisfies the usual flat spacetime Dirac equation. L Q ^ and k a are the angular and linear momentum 
operators of the particle. It follows from (II. 6 ) and (II.4) that the solution of (II.1) can be written in the form 

T(ar) = e -i$s {-i^(x)V M - to) e" i$G T 0 (a;), (II. 8 ) 

and also as 

^O) = - 7 T - - m ) e~ l ® T ^o{x ), (II.9) 

2 m 

where = <I> S + 4 >g is of first order in 'y a p(x). The factor —1/2 m on the r.h.s. of (II.9) is required by the condition 
that both sides of the equation agree when the gravitational field vanishes. 

It is useful to re-derive some known results from the covariant Dirac equation. On multiplying (II.1) on the left by 
{—i^ v [x)T> v — to) and using the relations 

v # ,r v (®) - v„r„(s) + i[r M (x), r„(®)] = ~a a 0 {x)R a ^ v , (ii.io) 
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and 


= - l -(J ap {x)R a ^ v , ( 11 . 11 ) 

we obtain the equation 

(g^V^ - | + m 2 ) (a) = 0. (II.12) 

In (II.11) and (11.12) Rap^u is the Riemann tensor, R the Ricci scalar, and <r a P(x) = (f/2)[7 Q (;r), 7 ^ (a)]. 

By using Eq. (II.4), we also find 

{-i^(x)D v - m) S (* 7 "V m - m)*(x) = S ( 5 ^V M V y + ?n 2 )4'( x) = 0. (11.13) 

While Eq. (11.13) is mainly a re-statement of the fact that (II. 8 ) is a solution of (II.1), Eq. (11.12) implies that the 
gyro-gravitational ratio of a massive Dirac particle is one when R ^ 0, as found by Oliveira and Tiomno [ 8 ], Audretsch 
[9] and Kannenberg [10]. 

It is known that the weak field approximation does not fix the reference frame completely. The 

transformations of coordinates x M —> x^ + with £ M ( x ) also small of first order, are still allowed and lead to the 

’’gauge” transformations 7 M „ —> It is therefore necessary to show that in (II.9) is gauge invariant. 

In fact, on applying Stokes theorem to a closed spacetime path C and using (II.10), we find that 4 >t changes by 

A$ t = \ [ dr^J^R^ap, (11.14) 

where E is a surface bound by C and J a @ is the total momentum of the particle. Eq. (11.14) shows that (II.9) is gauge 
invariant and confirms that, to first order in the gravitational field, the gyro-gravitational ratio of a Dirac particle 
is one even when R = 0. Use of (II.9) automatically insures that particle spin and angular momentum are treated 
simultaneously and correctly. 


III. NEUTRINO GEOMETRICAL OPTICS 

In this Section we study the propagation of a one-flavor Dirac neutrino in the Lense-Thirring gravitational field [11] 
represented, in its post-Newtonian form, by 


700 = 20 , 7 ij = 2(j)Sij , 70* = hi = -3 (J A r);, 


(III.l) 


where 


<t> = 


GM 

r 


h = 


4 GMR 2 lo , 

— 575 — 


~ x i 0 ), 


(IH-2) 


and M, R, w = (0,0, w) and J are mass, radius, angular velocity and angular momentum of the source. This metric 
is particularly interesting because it describes a rotating source and has no Newtonian counterpart. By using the 
freedom allowed by local Lorentz transformations, the vierbein field to 0( 7 ^) is 

e? = 0, e° = l — 0 , ei=hi, e[ = (1 + 0 ) 6 l k . (HL3) 

It is also useful to further isolate the gravitational contribution in (III.3) by writing ~ 0] + /i((. The spin connection 
can be calculated using (II.2) and (III.3) and is 


r 0 = -Io„ - '-l I,,, - (HI.4) 

r i = -\(Kj + h j,r)^ 3 - a' 3 ■ 


Explicit expressions for T M are given in Eq.(A.4) of the Appendix. 
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In the geometrical optics approximation, valid whenever \dij /JM \ <C kj IJM \. where k is the momentum of the particle, 
the interaction between the angular momentum of the source and the particle’s spin vanishes. This interaction is 
quantum mechanical in origin. Then the geometrical phase Tq is sufficient to reproduce the classical angle of deflection, 
as it should, because is common to the solutions of all wave equations when the spin is neglected. 

The deflection angle p is defined by 


X&iup 


~9ijP±Pl 


P II 


(III.5) 


where pj_ and p\\ are the orthogonal and parallel components of the momentum with respect to the initial direction 
of propagation. In the weak field approximation tan ip ~ <p and (III.5) reduces to 


ip ~ 


Ipj-I 

k\\ 


(III.6) 


where fcy = P\\ is the unperturbed momentum and |pj_| = y —'r)ijP l L p’ J _, for p l ± ~ 0( 7 ^,,,). 

It is clear from (II.7) and (II. 8 ) that, once To (a;) is chosen to be a plane wave solution of the flat spacetime Dirac 
equation, the geometrical phase of a neutrino of four-momentum k M is given by 

v(x) = —k a x a — $g{x) , (III.7) 

where 

$g(x) = -^J dz x [j a \^(z) - 7/?A,a0z)] ((x a - z a )kP - (x 0 - zP)k a ) + ^ J dz x j a xk a . (III. 8 ) 

The components of pj_ can be determined from the equation 

dv 

Pi = = -h - $G,i = (III.9) 

1 1 f x 

= ~h - - 7 ai(x)k a + - J dz x (-yi\, 0 (z) - 7/?a ,i( z ))k p ■ 

We consider below the two cases of propagation along the z-axis, which is parallel to the angular momentum of the 
source, and along the ai-axis, orthogonal to it. In both instances, the neutrinos are assumed to be ultrarelativistic, 
i.e. dz° ~ dz( 1 + m 2 /2E 2 ), E ~ k( 1 + m 2 /2E 2 ). 


A. Propagation along z 


Without loss of generality, we consider neutrinos starting from z = —oo with impact parameter b > R and 
propagating along x = b, y = 0. From (III.9), (III.l) and (III.2), we find 


P l = 


1 


= —2k 1 + 


[ dz°"/oo,ik° + j dz 3 733 ,ifc 3 

J— OO -J — oo 


(III. 10) 


2E 2 J J_ c 
1 
2 


4>,idz , 


! 1 r n 

P2 = -^702 k° + - d2°720,3fc 3 = 0 

J — oo 


and 


(p±y 

{PX? 


= g l ^Pv ^-pi = - 
= g 2li Pii — h 2 E = - 


2 GMk 




AGMR 2 ubk 

5?’ 3 


1 + 


m‘ 


2Ei) ' 


(III.ll) 
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From (III.6), we finally obtain 


2 GM 


ip = 


1 + 




2\2 f2R 2 b 2 0j\ 


b v ' 2E 2 ) \V ' r) ' V 5r 3 J ’ 


(III. 12) 


which is the deflection predicted by general relativity for photons, with corrections due to the neutrino mass and the 
rotation of the source. In the limit z —* oo Eq. (III. 12) reduces to 


p = 


AGM / 


1 + 


TO 


b V ' 2 E 2 ) ' 

If, as in [1], we introduce an effective mechanical momentum 

Pa f = k a + ®G,a + r Q , 


(III. 13) 


(III. 14) 


then there is a residual contribution of the real part of F 2 to ip. It is quantum mechanical, it amounts to 
' GM 5 % ~ b ~ ^1 — and is smaller than all the other terms in p. With this correction, we obtain 


p = 


2 GM 


m 

E 2 




5?’ 3 


5 Er 5 


1 - 


E A 


(III. 15) 


B. Propagation along x 

In this case neutrinos start from x = — oo with impact parameter b. For simplicity, we consider neutrinos propagating 
in the equatorial plane z = 0, y = b. We find 

P 2 = -^ 702 ^° + \ J dz° ( 720 ,lfc 1 - 700 , 2 ^° - 710 . 2 & 1 ) - ^ J dz 1 ( 710 , 2 k 0 + Jupk 1 ) (III.16) 

= ~ k (y 1 + ^ 2 )/ ( 2 ^,2 + h lt2 )dx , 


P3 = 0 


(III. 17) 


and 


(p±) 2 = g 2tl Pii - -P2 + h 2 E = - 


2 GMk 


1 + 


to 2 \ / 2R 2 u 

2E 2 ) V 5 V 




(p±) 3 = g 3fJ, p l i - -P3 = 0 . 

It then follows that 

2 GM 


( 2R 2 lu\ ( m 2 \ ( x\ 

{ —) ( 1 + 2 &) ( 1 + 7 ) • 


(III. 18) 


(III. 19) 


Contrary to the previous case, the contribution of the angular momentum of the source does not vanish in the limit 
x —> 00 . In fact, in this limit we get 


p = 


AGM 

b 




(III.20) 


which coincides with the prediction of general relativity. Additional, smaller spin contributions can be obtained from 
as in Subsection (III.A). 
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IV. HELICITY TRANSITIONS 


In what follows, it is convenient to write the left and right neutrino wave functions in the form 


’I'o(z) = J'OL.fle 


—ik a x 


E 


m 


2 E 


vl,r 

<t k 

E+m L,R 


„—ik n 


(IV.l) 


where cr = (a 1 , a 2 , cr 3 ) represents the Pauli matrices. vl, R are eigenvectors of er • k corresponding to negative and 
positive helicity and Dq L,Ft{k) = Vq l R (k) 7 ®, R {k)^o L,n(k) = 1. This notation already takes into account the fact 

that if v± are the helicity states, then we have ~ is-, v R ~ v + for relativistic neutrinos. 

In general, the spin precesses during the motion of the neutrino. This can be seen, for instance, from the contribution 
<& s in The expectation value of the contribution of To to the effective mechanical momentum (III.14) can in fact 
be re-written in the form 

(IV.2) 


where fl = G ^ r 3 ^1 — j cD. Eq.(IV.2) represents the spin-rotation coupling, or Mashhoon term [12], for the 
Lense-Thirring metric. Here rotation is provided by the gravitational source, rather than by the particles themselves. 

We now study the helicity flip of one flavor neutrinos as they propagate in the gravitational field produced by a 
rotating mass. The neutrino state vector can be written as 


= a(X)\u R ) + P(X)\u L ), (IV.3) 

where \a\ 2 + \(3 \ 2 = 1 and A is an affine parameter along the world-line. In order to determine a and /3 , we can write 
Eq. (II. 8 ) as 


where 


MA)> = f(A)hMA)), 


T = -E (-^(xyDp - m) e I<t>T , 


(IV.4) 
(IV.5) 
(IV. 6 ) 


and |^o(A)) is the corresponding solution in Minkowski spacetime. The latter can be written as 

IV’o(A)) = e~ lk ' x [a(0)|^ fl ) + P(0)\i/ L )] . 

Strictly speaking, |^>(A)) should also be normalized. However, it is shown below that c*(A) is already of 0( 7 M „), can 
only produce higher order terms and is therefore unnecessary in this calculation. From (IV.3), (IV.4) and (IV. 6 ) we 
obtain 


a(A) = (z'kIV’(A)) = ot(0)(u R \f\u R ) + P(0){v R \T\v L ). 


(IV-7) 


An equation for /3 can be derived in an entirely similar way. 

If we consider neutrinos which are created in the left-handed state, then |a(0 )| 2 = 0, |/3(0 )| 2 = 1, and we obtain 


Pl->r = |«(A)| 2 = {v r \T\vl) 


' Ao 




(IV. 8 ) 


where = k^/m. As remarked in [ 1 ], x M need not be a null vector if we assume that the neutrino moves along an 
’’average” trajectory. We also find, to lowest order, 


d^f = 2 - ^-i2?n4> GjM - i{j a k a + m)$ s ^ + 7 + $G,a/i)) 
$ s ,a = T a , 4>G,a M = fc^Tf , ^o(7 Q fc Q +m) = 2E^7 6 , 


(IV.9) 


where are the usual Christoffel symbols, and 


(v R \x ,i d IJ f\p L ) = 


E 




(IV. 10) 


In what follows, we compute the probability amplitude (IV.10) for neutrinos propagating along the z and the x 
directions explicitly. 


6 




A. Propagation along z 


For propagation along the 2 -axis, we have k° = E and k 3 = k ~ E( 1 — m 2 /2E 2 ). As in Section III, we choose 
y = 0, x = b. With the help of Eq.(A.7), we get 


rv rv^ 

2 mE 


,k x _ k m 

-*— v r T x v l = —0,i + *—«2,3 , 
m m 4 b 

(+> + +)47 4 ^ = -A(i + V) . 


(IV.ll) 


Summing up, and neglecting terms of 0(m/E) 2 , Eq.(IV.lO) becomes 

{u R \x^d^f\u L ) = + l -h 2 , 3 ■ 

The contributions to 0((E/m) 2 ) vanish. As a consequence 

da m da m (1 i \ 

~dz ~ ~E~d\ = ~E (+•' + W 3 ) ' 

and the probability amplitude for the vl transition is of 0(m/E), as expected. 

Integrating (IV.13) from —oo to z, yields 


(IV. 12) 


(IV. 13) 


E 


m GM 

~E~2b~ 


I dz</) t i + jh 2 {z) 

^ 2 ujR 2 b 2 ' 

1 r ^ * br 3 


(IV. 14) 


It also follows that 


Pl^r{- 


fm\ 2 (GM\ 2 [/ 2;\ 

hr ( 1 + t) 


2 ub 2 R 2 
5 r 3 


(IV.15) 


The first of the two terms in (IV. 15) comes from the mass of the gravitational source. The second from the source’s 
angular momentum and vanishes for r —» oo because the contribution from —oo to 0 exactly cancels that from 0 to 
+oo. In fact, if we consider neutrinos propagating from 0 to +oo, we obtain 


Pl- 


(m \ 2 ( GM\ 

^He) hr) 


i+ 


(2 ujR 2 
V 5 b 


(IV.16) 


According to semiclassical spin precession equations [13], there should be no spin motion when spin and uj are parallel 
as in the present case. This is a hint that rotation of the source, rather than of the particles, should produce a similar 
effect. The probabilities (IV.15) and (IV.16) mark therefore a departure from expected results. They are however 
small of second order. Both expressions vanish for m —> 0, as it should for a stationary metric. In this case, in fact, 
helicity is conserved [14], It is interesting to observe that spin precession also occurs when uj vanishes [15,16]. In 

the case of (IV.15) the mass contribution is larger when b < { r /R)\J~^;■> which, close to the source, with b ~ r ~ R 1 
becomes Ruj < 5/2 and is always satisfied. In the case described by (IV.16), the rotational contribution is larger if 
b/R < 2u>R/5 which effectively restricts the region of dominance to a narrow strip about the 2 -axis in the equatorial 
plane, if the source is compact and u> is relatively large. 


B. Propagation along x 

In this case, we put k° = E, k 1 = k ~ E( 1 — m 2 /2E 2 ). As in Section III, the calculation can be simplified by 
assuming that the motion is in the equatorial plane with 2 = 0, y — b. We then have 
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h X h 

A/ rv^i 

2m E 


k x _ ^ k E 2 -\-k 2 E 2 — k 2 

~l VrY\Vl = l — (p, 2 + l— 7 —^ 1,2 % — —h 2,1 

m m AmE AmE 

«, + (i + i) - f 2si ft w • 


Summing up, and neglecting terms of 0(m/E) 2 , Eq.(IV.lO) becomes 

{iy R \x^d^T\iy L ) = ~0 )2 + j(/ii,2 - h 2 ,i) ■ 
The contributions to 0((E/m) 2 ) again vanish and we get 


da m da m 
dx E dX E 

Integrating (IV.19) from — oo to x, we obtain 


X % 

^,2 + 4 (^ 1,2 - ^ 2 ,l) 


0(m/E). 


. m GM 


1 - 


2 ojW 


5 b 




and 


Pl->r{~ 00 , x) ~ 




GMY 

\Tm) 


1 - 


2u)R* 


5 b 


( 1 + - 


) 2 

r / 


(IV. 17) 


(IV. 18) 


(IV. 19) 


(IV.20) 


(IV.21) 


Obviously, the mass contribution is the same as for propagation along the z-axis. However, the two cases differ 
substantially in the behavior of the angular momentum term. In this case, in fact, this term is even, so it does not 
vanish for r —» 00 . If we consider neutrinos generated at x = 0 and propagating to x = + 00 , we find 


P^O.+oo) =/(=/( 


GM 


1 - 


2 wl? 2 
5 b 


(IV.22) 


The mass term is larger when 2u i^ <1. At the poles b ~ R and the mass term dominates because the condition 
uR < 5/2 is always satisfied. The angular momentum contribution prevails in proximity of the equatorial plane. The 
transition probability vanishes at b = 2uR 2 /5. 


V. FLAVOR OSCILLATIONS 


Eq. (II.9) can be recast into a form that contains only first order contributions. By using (II.9) we find 

*(*)=/( aOe-^- ife “*“To, 


(v.i) 


where 

f{x) = 2- [e^ 7 & (fc M + $g, m ) +m] (V.2) 

and H/o now represents the phase independent part of (IV.1) that refers to mass eigenstate neutrinos. The relationship 
between flavor eigenstates (Greek indices) and mass eigenstates (Latin indices) is given by the standard expression 

KW) = E,U aj {0)\v j (x )}, (V.3) 

3 


into which (V.I) must now be substituted. We find 


K0r(A))> = 


Y,U aj O(x( A))K), 

3 


(V.4) 



where 


O(A) = ^ Kl & {K + *g,m) + m] e-^ T+k ^ , 

and U is the mixing matrix 

cos 9 sin 9 
— sin 9 cos 9 




(V.5) 


(V.6) 


Restricting a to the flavor e for simplicity, we define the column matrix 


X = 


= C/t 


(l/ e \Ue(x)) \ _ 
{Vy\v e {x)) 

On Oii 

021 022 


cos 2 80n + sin 9 cos 90 12 + sin 9 cos 0 O 2 i + sin 2 00 22 
— sin 9 cos 90n — sin 2 90i2 + cos 2 f?0 2 1 + sin 9 cos 0O 22 


u<!>. 


where Oy = (^|0(A)|^) = = (Ci|0(A)|Cj)^ij, i,j = 1,2, and </> 

The differentiation with respect to the parameter A gives 


1 

0 



c/t 


6 dA 


0 


• dp 2 2 
6 dA 




(V.7) 

(V.8) 


(V.9) 


By keeping only terms of the first order in j a p, observing that $o e a P ~ 0(7#L)> that 4 >g, ^r p « 0( 7 2 „) and using 
(A.2) we find that the matrix elements idOij/dX in (V.9) are of the form 


. dOij 

1 - - 

dX 


A« Oij + 'EdflOu+C, 


(i) 


(V.10) 


where the index (i) refers to the i-th mass eigenstate, and 


A® = x (i)p 


Ox 


+ $ 


« 

G,p 


= (v i \x p T p \v j )=xM p {QT p \Q)5 i] , 


(V.ll) 




.X 


(i)p 




2 rrii 


KpW 




„ — —ik^ • x 

; e 0 


The equation of evolution therefore is 


where 


and 


l -± = (A f +df)x + C, 

A <= u '{ T ) u * 


c = u f 


0n 0 
0 0 22 


U<t>. 


(V.12) 


(V.13) 


(V.14) 


A relation similar to (V.13) links df to d. 

Up to (V.14) the derivation has been general and can be applied to any metric. It is not possible to say more about 
the matrix 0 at this point, except that it may contain some dissipation terms. 

In the case of the Lense-Thirring metric, however, 0 can be written as C = C\- After some algebra, the matrix 
elements 0y become, in fact, 

C ’ij ~ -*f'7oo,3e- i(fc(i> ' (x - xo) -*o )) <5« - (V.15) 
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Starting from Eq. (V.5), we can write 


Oij ~ e -ifc (i) -(x-*o (1 + A (i) + i2J 0 )S i:j , 


(V.16) 


where 


A + = \K 


(V.17) 


Jn = 


GMR 2 u 1 
10 x 2 + y 2 


/£\ 3 _ f £o V 
\rJ \ r o J 


GM xy 


4 r \x 2 + z 2 y 2 + z‘ 


To O^oom 2 /E 2 ), = A is independent of neutrino mass and energy. It is, in fact a pure geometrical term. Using 

(V.16), Eq.(V.15) becomes 


where 


ri\ % ) _ r< 

°ij ~ U U 


■ E Too,3 
4: 1 + A + z2 Jq 


Oij = COij 


c 


,E 

-1-^ 700,3 • 


To 0{m 2 /E 2 ), C does not depend on the mass eigenstate index (i). Eq. (V.10) therefore becomes 


and the matrix C is 


. dOij 

i - - 

d.X 


A® Oij+Y,4kO kj +COij, 

k 


C=U t 


CO 11 0 

0 CO 22 


U(j>= Cx • 


(V.18) 


(V.19) 


(V.20) 


(V.21) 


In the equation of evolution (V.12) one can then carry out the transformation y —> e 1 ^ x o Cd> ' y and remove the C 
dependence by means of a suitable normalization. 

The equation of evolution can be re-cast into a more traditional form by taking x p as the tangent vector to the null 
world line. Then x p x p = 0, k p = (E, k l ) ~ (E, n l E(l — m 2 /2E 2 )), where n l is a unit vector parallel to the neutrino 
three-momentum, k° = x° and k l ~ (1 — e)x % , with £<1. The matrix elements of (V.13) contain terms of the form 

~ _|_ sE 2 . The term with E is diagonal in the matrix of evolution and does not therefore contribute to 

the oscillations. The equation of evolution (V.9) then becomes 

= {jY + ^G+^Xi (v.22) 


where the flavor mass matrix Mf is related to the vacuum mass matrix in the flavor base by 


while 


and 


1 0 (Czlk^aM ) L ' 


T f = U f ^ (Cil^pICi) 


(C2 1 fc p T p | C2) 


U. 


(V.23) 


(V.24) 


(V.25) 
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To 0('i liV m 2 /E 2 ), Eq. (V.22) becomes 


- (^2 + kP{ ^» + r ^ )} ) X • (V ' 26) 

The term in brackets on the r.h.s. of (V.26) is the effective Hamiltonian of evolution. 

By writing the effective mechanical momentum in the form (III. 14), squaring it and keeping only first order terms, 
we find 


and 


M 2 ~ m 2 + 2/ c q $ t , « 


(V.27) 


M 2 


\peff\ = ^ E 2 _ M 2 ~ E - =E - 

1 1 V 2 EE 


($G, a +T a ) 


(V.28) 


The third and fourth terms of (V.28) represent the gravitational contributions (V.24) and (V.25) and provide an 
explicit expression for the term p ■ AfV^ in Eq.(27) of [1] when matter effects are neglected. These can be easily 
incorporated in the diagonal element of (V.22). The effective potential V v . induced by matter, depends on neutrino 
flavors and is defined by [17] 


K. = -V Pe = V 0 (3Y e -1 + 4 Y„J, 
K MlT =-V^ r = V 0 (Y e -l + 2Y„ e ) f 


where Y e (Y v ) represents the ratio between the number density of electrons (neutrinos), 

Gpp p 


Vo = 


3.8eV, 


V2m n 10 14 gr/ cm 3 

m n = 938 MeV is the nucleon mass, p the matter density, and Gf the Fermi coupling constant. 


(V.29) 

(V.30) 


(V.31) 


A. Transition amplitudes 

In order to calculate the transition amplitudes between states of different flavor, we follow a procedure analogous to 
that developed in Section IV for the helicity transitions. To discuss oscillations, we assume that neutrinos in different 
mass eigenstates have the same energy and keep in mind that energy is conserved in a stationary gravitational field. 
In order that different neutrinos interfere at the same final point, we require, as usual, that the relevant components 
of the wave function do not start from the same initial point. This however introduces negligible corrections in the 
final results if one assumes, as we do, that Am 2 = m 2 — m 2 <C m 2 and 7 ^ <C 1. 

Let us consider, for simplicity, only two flavors and, in particular, the transition u e —► v^. The transition probability 
is given by 


P,, 


<^ MA ))| 2 


MA)k(A)) ’ 

where |z/ e (A)) represents at A a neutrino initially created in the state \v e ). From (V.3) and (V. 6 ) we obtain 

(^lz'e(A)} = sin 0 cos 6 >((i/ 2 |^ 2 (A)) - (j'ilz'i(A))) 

(i/ e (A)|i/ e (A)) = cos 2 6{v\(X)\v\(X)) +sin 2 0 (j/ 2 (A)|^ 2 (A)). 

For a generic mass eigenstate, we find 

^Hz'(A)} = ^ = —ix ,i kn(u\u(X)) + x* 1 {v Q \d ll f\v 0 )e-' lk > 1 - x> '. 

Since {vQ\d^T\v 0 ) is already ~ 0( 7 M „), it is possible to carry out the substitution 


(V.32) 


(V.33) 


(V.34) 


e-**** = {v\MX))* WWW)- 


(V.35) 
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Eq.(V.34) has the first order solution 


WWW) - exp 




where P and Q are the points at which neutrinos are generated and detected. 
We also obtain 


= x» (voider+ T')Wo), 


(V.36) 


(V.37) 


which has the first order solution 


From (IV.9) we get 


(z'(A)l^(A)} ~ exp / J d\x^Wld^T + T t )|j/ 0 )j • 


Wo\dfj,TWo) — + r^J 

(u 0 \d^f + f^Wo) = ^(h^ + r^). 


Substituting (V.39) into (V.32), we finally obtain 


Pu. 


( [ Q E _ 

sin 9 cos 9 exp — i / dXx ^— 

V Jp rn 


exp 


~i J p d\x»(kW + - ex P J p dXx^ik^ + $^)j 

= sin 2 26 sin 2 J dXx^(Ak fJj + . 


(V.38) 


(V.39) 


(V.40) 


In deriving (V.40) we have used the fact that the terms i7 0 a^ u Po are always proportional to m/E, as shown in the 
Appendix, and that the contributions from the coefficients F M are therefore independent of to, E and k. It follows 
from (V.40), that the effect of the gravitational field is simply obtained by the substitution k^ —> /c p + $g> — IV 
This conforms to the semiclassical principle that action is the quantum phase of the particle. 

Eq.(V.40) can be simply rewritten as 


P Ve ^ Vil = sin 2 29 sin 2 



dx-’ (Akj + 


(V.41) 


because the neutrinos have fixed energy, hence Afco = 0, and <I>g,o = 0 for static gravitational fields. 

If, for simplicity, we consider the motion to be parallel to the aA-axis, the first term can be written, to 0(m/E), in 
the traditional form 


1 Am 2 

2 Jp o). 


Using (III.8), (III.9) and dz° = ( E/k)dz l , we can write 

1 1 f z 

®G,i( z ) = “ 2 J dzX ^ iX1/3 ~~ 7/3A ’0 


TO 


— E(2</>(z) + hi(z)) + + E ( 1 + 


TO 

2E* 


<f>( z o), 


(V.42) 


(V.43) 
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where, as above, we have retained only the lowest non vanishing terms in m/E. We finally obtain 


A* O .<W = ^W + ^(--0>- (V.44) 

The first term of (V.44) is the well known gravitational contribution to the phase of oscillating neutrinos. The second 
term corresponds to the redefinition of the constant E —» E^Jgoo{zo) in the non-gravitational term of the oscillation 
for a neutrino created at zo, at the desired order of approximation. With this redefinition, the phase that appears in 
(V.40) becomes 


n = 


Am 2 
4 E 2 


(z - z 0 ) + 


m 2 Am 2 
2 E 3 



(V.45) 


The coordinate difference z — Zq is not, however, the physical distance between two points. This is in fact defined as 
l = \J —gadx % ~ J“ (1 — (f>)dx l : while E is the neutrino energy measured by an inertial observer at rest at infinity. 

Introducing Ei = Egoo(z), which is the energy measured by a locally inertial observer momentarily at rest in the 
gravitational field, the phase O can be entirely rewritten in terms of physical quantities as 


n = ^ [ l ^L 

4 Jo Ei 


to 2 Am 2 
2E 3 



(V.46) 


Eq.(V.46) reflects the fact that the curvature of spacetime affects the oscillation probability through the gravitational 
red-shift of the local energy Ei and the proper distance dl. 

The oscillation amplitude does not depend on 70 *. This means that the relevant oscillation parameters like the 
phase fi, depend quadratically (and not linearly) on the angular momentum of the source. This should be expected 
because, as noted by Wudka [18], the quantum mechanical phase is a scalar, whereas the angular momentum is a 
pseudovector. 


VI. CONCLUSIONS 

In this paper we have applied the weak held solution of [2] to the study of neutrinos in a Lense-Thirring Held. 
There are advantages to using this solution. It is exact to first order, is covariant and gauge invariant. It is based on 
the Dirac equation and reproduces well all gravitational effects that have so far been observed [19,20], or have been 
predicted to exist for a spin-1/2 particle [6,7]. It also refers to fermions with unit gyro-gravitational ratio [3] and 
thus allows a unihed treatment of spin-rotation and angular momentum-rotation coupling without requiring ad hoc 
procedures. 

We have applied the solution to the propagation of single flavor neutrinos. Two cases have been considered in the 
geometrical optics approximation. In both instances the larger contribution to the deviation ip comes from the mass of 
the source and is as predicted by Einstein’s theory, with corrections due to the neutrino mass. For propagation parallel 
to the axis of rotation of the source, the rotation corrections vanish at infinity. Not so for propagation perpendicular 
to the axis of rotation. 

It may be argued that spin contributions exist in the approximation used by introducing an effective neutrino 
mechanical momentum given by (III. 14) that incorporates a contribution from T^ as in [1]. This definition is also 
supported by the form of the flavor oscillation Hamiltonian on the r.h.s. of (V.26)-(V.27). The contributions that 
come from taking the expectation values of T^ are given in (III.15) and are much smaller than those of (III.12). 

We have then calculated the helicity transition amplitudes of ultrarelativistic, single flavor neutrinos as they prop¬ 
agate in the Lense-Thirring Held. These transitions are interesting because at high energies chirality states are 
predominantly helicity states and right-handed neutrinos do not interact. The transition probabilities are of 0( 7 2 „). 
Two directions of propagation have again been selected and the results contain contributions from both mass and 
angular momentum of the source. The transitions also occur in the absence of rotation or with spin parallel to 
rotation, which is unexpected on semiclassical grounds. The mass contributions predominate when the neutrinos 
propagate from r = 0 to r = 00 (and matter effects are neglected), provided b > 2wi? 2 /5. There is, however, a 
narrow region about the axis of propagation in the equatorial plane where the u> contribution is larger. The rotational 
contribution behaves diHerently in the two cases. It vanishes as z —» Too for propagation along 2 , but not so as 
x —> 00 in the second case. In addition, when the neutrinos propagate from x = 0 to x = 00 , the mass term dominates 
in the neighborhood of the poles, while the contribution of u> is larger close to the equator, with no attenuation at 
b = 2 ujR 2 /5. 
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We have also calculated gravity induced, two-flavor oscillations and derived the relative equation and effective 
Hamiltonian. A comparison with the simplified approach of [1] leads to explicit expressions for the term p ■ AfVc of 
these authors. The expressions contain F^, as expected, and also These terms do not exhaust, however, all 

possibilities if the metric is not stationary, because of the presence of the matrix C in (V.12). In the general case, 
therefore, the evolution of \ is more complicated and the effective Hamiltonian does not comply with the simpler form 
given in [1], or in (V.26). Finally, the transition probabilities do indeed oscillate for the Lense-Tlrirring metric, and 
the curvature of spacetime enters the oscillation probability through the gravitational red-shift of the local energy Ei 
and the proper distance dl. 

The neutrino mass appears quadratically in all effects considered and, in particular, as Am 2 in flavor oscillations. It 
therefore seems possible, in principle, to obtain the neutrino masses from an appropriate combination of measurements 
of these effects, without resorting to models that might depend on additional aspects of particle physics. The only 
ingredient used is in fact the covariant Dirac equation whose merits have been extolled at the beginning of this section. 

The results presented in this paper agree with those of other authors, where appropriate. They can be applied to 
a number of problems in astroparticle physics and cosmology [21]. For instance, an interesting question is whether 
gravity induced helicity and flavor transitions could effect changes in the ratio v e : : v T of the expected fluxes at 

Earth. 

Lepton asymmetry in the Universe [21] also is an interesting problem. It is known that the active-sterile oscillation 
of neutrinos can generate a discrepancy in the neutrino and antineutrino number densities. The lepton number of a 
neutrino of flavor / is defined by Lf = (n V} — n P/ )/n 7 (T), where n Uf (n^ f ) is the number density of neutrinos (an¬ 
tineutrinos) and ?r 7 (T) is the number density of photons at temperature T. As discussed in Sec. IV, the gravitational 
field generates transitions from left-handed (active) neutrinos to right-handed (sterile) neutrinos. If, in primordial 
conditions, (IV.16) and (IV.22) become larger, then helicity transitions may contribute in some measure to lepton 
asymmetry. 


APPENDIX A: USEFUL FORMULAE 


In this Appendix, we collect useful formulae and results that have been extensively used in the paper. 
Derivatives of $< 3 . The first derivative with respect to gives 

dz x (j l _ lX ,0 - lp\^)k 0 + ^ 7 a^k a , (A.l) 

whereas the second derivative is 




flV 5 


(A.2) 


where rare the Christoffell symbols of the second type. 

Christoffel symbols and spin connections. To 0( 7 ^), the Christoffel symbols for a Lense-Thirring metric are 


r oo = 0, = <j > ti , T% = l -{h^ + h hi ), 

r oo = ^U r o j = ~ bid) > = 

The spin connection coefficients, already calculated in (III.4), have the explicit form 


r n = - 


GM 


r, = 


r 2 = 


r, = 


2 r 3 

3 GMR 2 u r. 

5r 5 
3 GMR 2 u r 
5r 5 

3 GMR 2 lo 
5 r 5 


xa 01 + ya 02 + za 03 ) + 


gmr 2 lu r 


5 r 5 


t 03 


2 xya m + (y 2 - ar)er u “ + yza 1 

( y 2 - x 2 )a - 2 xya 32 - xza 33 
GM 


(r 2 - 3 *> 12 + 3 yza 13 - 3 xza 23 
GM 


yza 01 - xza u * ) + 




xa 13 + ya 


2 r 3 
GM 
2?’ 3 
23 ^ 


(ya i2 + za i3 ) 


'^xa 12 + za 23 ' 


(A.3) 


(A.4) 
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7 and a matrices. We have used the Dirac representation. The neutrino mass eigenstates are 


v(x) = v 0L e 


—ik a 


E + m 


vl,r 


erk 


—ika 


(A.5) 


2 E V e+Ei u l,r 

Motion along z direction. In this case vl, = ^ ^ ^ and Ur = ^ q ^ (he. v l,r are eigenstates of a 1 2 3 4 5 ). To calculate 
the helicity transitions, we need 


{"Rh °v 01 W L ) = ~i^, {vR\'f<J 02 \v L ) = -^, {vR.\y Q cr m \v L ) = 0 , 

{vr\^(j V2 \u l ) = o, (RRh°<j 13 W L ) = i, Wl7 6 cr 23 K) = 1 , 

{vr\7°Wl) = Wr\i 3 Wl) = 0 , (rrI^Wl) =, {vrI^Wl} = ■ 


(A. 6 ) 


As a consequence of (A. 6 ) and (III.3) we obtain 

.,6 


k 1 

(z'fl|7 U ro|z'L) = ^(i^, 1 + 0,2) — — ^3,1) + (/i2,3 — ^3,2)] 

k 1 

(Fr|7 TsI^l) = ^[i(/ll,3 + ^3,1) + (/l2,3 + ^3,2)] + ^(*0,1 + 0,2) ■ 


(A.7) 


find 


Motion along x direction. In this case 17 , = -i= ^ ^ ^ and jzr = ^ j ^ (i.e. vl.r are eigenstates of a 1 ). We 

{vr\i°o- 01 \v l ) = 0, (^r|7 0 ct 02 |i/z,) = , (A.8) 

(l/ fl |7°cr 12 |^ L ) = 1, (vr^o-^Ivl) = -i, (vr\^(J 23 \v L ) = 0, 

(z'aIt'Vl) = (^rItVl) = 0, {vrY^Wl) = , Wl7 3 ki) =--| • 


We also obtain (see (III.3)) 


(*7?l7° r oK) = ^(-<^,2 + #,3) ^ 7 [(hi, 2 - h 2 , 1) - *(fti,3 - *3,1)] 


2# 

k 


(A.9) 


(^k| 7 r i|^i) = Jg[-(/ll,2 + /l2,l) +i(ftl,3 + * 3 , 1 )] + o(—0,2 +#, 3 ) ■ 
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